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Abstract: We show for a large class of interacting particle systems that whenever the 
stationary measure is not reversible for the dynamics, then the mean entropy produc- 
tion in the steady state is strictly positive. This extends to the thermodynamic limit 
the equivalence between microscopic reversibility and zero mean entropy production: 
time-reversal invariance cannot be spontaneously broken. 

Keywords: stochastic interacting particle systems, entropy production, (generalized) 
detailed balance. 



1 Introduction 



Reversibility and entropy are words with many meanings even within the context 
of nonequihbrium statistical mechanics. One class of models that has often been 
considered for learning about nonequihbrium behavior is that of interacting particle 
systems. These are stochastic dynamics for spatially extended systems in which 
particles locally interact. They are mostly toy- models remaining far from realistic in 

their microscopic details. Yet, it is believed that for some good purposes, the details 
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do not matter so much and one should be concerned more with the symmetries, 
possible conservation laws, locality of the interaction etc. to hope to understand 
something about real nature. 

This paper is about the relation between time-reversal invariance and the positivity of 
entropy production. We do this in the context of interacting particle systems following 
the work in [|10|, ||, |^, ^ 0|. The physics background will be discussed in Section 



3. The main question is to understand why there cannot exist a 'superconducting' 
interacting particle system in the sense of the title of this paper to be specified below. 
To understand the mathematical problem, let us look first at a finite Markov chain. 
Suppose that K is a finite set on which we have an involution tt : K ^ Kjir"^ = id, 
called time-reversal. Often, the most natural choice for interacting particle systems 
is TT = id because we think of the state space as consisting of occupation variables or 
of classical spins but our mathematical set-up will be more general. 
Let {Xt, t G [—T, T]) be a stationary Markov process (steady state) on K with law Pp. 
The subscript refers to the unique stationary probability measure p on we assume 
that p(a) > 0, a G -ft'. The rate to go from a to 6 is denoted by k{a, b),a,b ^ K and 
we assume that k{TTb, vra) = iff k{a, b) = (dynamic reversibility). The generator is 



Lf{a) = J2kia,b)[fib)-fia)] (1.1) 
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The time-reversed process of (Xf) is the stationary Markov process {Yt,t G [— T, T]) 
on K with Yt = irX^t having transition rates 

~k{a,b) = k{7fb,7fa)($^ (1.2) 
p{TTa) 

We denote its law by Pp^ (pvr is stationary for (Yt)). Of course, it easily happens 
that p = pn and yet, Pp 7^ Fp^,. The corresponding generator for the time-reversed 
process is L = 7iL*7i where the * refers to the adjoint with respect to the stationary 
measure p. 

We say that the process (Xj)^j. is 7r-reversible if Pp = Pp^. This implies that 
the stationary measure p satisfies p = pir and 

p{a)k{a, h) = k{nh, ■Ka)p{h), a,b E K (1.3) 

which is generalized (or extended) detailed balance (microscopic reversibility). For 
the generators, we then have L = L. Observe that (|1.3| ) by itself implies that p{a) = 
pT{{a)p{h) I pTi{h) whenever k{a, h) 7^ 0. Applying this successively with 61, . . . ,bn E K 
for which k{a,bi), k{bi,b2), . . . ^k{bn.,T^a) 7^ 0, we find that p(a) = p7T{a). On the 
other hand, vr-reversibility implies that p = pn is stationary. 

The entropy production is the random variable obtained from taking the relative 
action on pathspace with respect to time-reversal, see [Q for a recent review. Let p-T 
be a probability measure on K which we use to sample the initial data at time —T 
for the stochastic time-evolution generated by L. The law of this process is denoted 
by Pp_j,. Suppose now that for this process the state at time T is described by the 
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probability measure pt- We could as well start our process (at time —T) from pt'k 
and then obtain the process Pp^^Tr- For a particular realization uj = {uj{t),t e [— T, T]) 
of this process we let O^o; = (7ra;(— G [— T, T]) be its time-reversal. The entropy 
production Rt^{L, p^T, T) is a function of the reahzation over the time- interval [— T, T] 
and is then obtained as 



Here we are only interested in its steady state expectation value, that is the mean 
entropy production rate, which in fact can be written as 



where lEp denotes expectation with respect to Pp. The notation MEPt,{L, p) reminds 
us that this number depends on the transformation tt, the dynamics (generated via 
L) and the stationary measure p. The mean entropy production thus measures the 
degree to which Pp can be distinguished from PpTp. The main property of the mean 
entropy production is then: 

Proposition 1: Consider the stationary process {Xt) above with p — pir. Then, 
MEPjr(I/, p) = MEP^(L, p) > with equality if and only if the process {Xt) is ir- 
reversible. 

This says that for finite state space Markov chains there can be no current without 
heat, meaning that detailed balance is equivalent with zero mean entropy production. 
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(1.4) 





The problem we address here is whether the same remains true in the thermodynamic 
hmit, that is for spatially extended interacting particle systems. In this case we really 
should be speaking about the mean entropy production density, i.e., per unit volume, 
but we will not use this extension. Note that in this case and from now on we will 
not and we cannot assume in general that p = pn even if both are stationary. 
We discuss the general physics set-up and further interpretations in Section 3, after 
stating our mathematical results in Section 2. We start however with three examples 
illustrating some aspects. 

1.1 Examples 

Example A: We consider particles hopping on the one-dimensional lattice with a 
preferred direction that is itself subject to independent flips. The state space is 
{-l,+l}x{0, Ij^and the process is determined by choosing a constant rate c{E, tj) = 
1 for changes from a configuration (-B, 77) to {—E,rf) and taking rates 

c{x, E, rj) = e^r]^{l - r?^+i) + e~^r]^+i{l - r]^) 

for changes to {E,r)^'^~^^) where {r)^'^~^^)y — r]y ii x ^ y ^ x + 1, and {r]^'^^^)y = r]x 
when y — x + 1 and = rj^+i when y — x. The resulting Markov process has generator 

X 

+ f{-E,rj)-f{E,7j) (1.6) 
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For invariant measure p we take 

p{E,dT]) = ^iSE,+i + Se,-i) X Uu{dT]) 

where is the Bernoulh measure with specified density u G (0, 1). For time-reversal 
we take 7T{E,r]) = {—E^rf) so that p = pvr. 

It is easy to see that the process satisfies generahzed detailed balance, like in 
the sense that p{E, drj) = p{—E, dif) = p{E, drj^'^^^) and both 

c{E, 7]) = c{-E, 7]) and c{x, E, r]) = c(x, -E, 

The last identity depends of course crucially on the fact that n is not the identity and 
reverses left and right as preferred direction. At the same time, as can be computed 
explicitly, the mean entropy production is zero. The same remains true for tt a 
particle- hole transformation, {7Tr])x = ^ — Vx, leaving the field E unchanged. Then, 
p ^ pn for u ^ 1/2 but still generalized detailed balance holds. Finally if, instead, 
we were to take vr = identity as time-reversal, then we break the detailed balance 
condition and we obtain a strictly positive mean entropy production. 
Example B: We take the simplest example of a spinfiip dynamics for which the one- 
dimensional Ising model is stationary but not reversible (for tt =id). Exactly the 
same can be done in two dimensions, see 0. Spinfiips are transformations Ux '■ cr ^ 
Ux{(y) = a^, X E 'E,aE {+1, —1} ^ for equal to a except at the site x. 
Consider the one-dimensional spinfiip dynamics with the following asymmetric rates: 

c(x, a) = exp{-2pax(Tx+i) (1.7) 
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The invariant measure p is the one-dimensional Ising model at inverse temperature j3. 
The process starting from p is not time-reversal invariant and the entropy production 
is equal to MEP(L,p) = 4/9tanh/3 (that is with time-reversal tt =identity). On the 
other hand, this time-reversed process is easy to find; it is a spinfiip process with 
generator 



Let us now take for time-reversal tt the refiection: (vrcr)^; = (T-x which leaves p 
invariant. Since 



vrcr = 7r o" 



L* = ttLti and we have in fact generalized detailed balance ( |1.3| ): 

{a)=e 



c{—x, (ttct) ^) dp 

The denominator in the left hand side is the rate in the original process by which 
ttUxO- = 7r(cr^) = (vrcr)""'' is changed to na. As a result, MEP^(L,p) =MEP.„{L*, p) = 
0. 

Example C: Instead of driving the system in the bulk and breaking detailed balance 
via some external fields that act on each component of the system, we may also 
consider boundary driven processes. For this we need to start with finite volume. 
The simplest interesting case is that of a symmetric exclusion process on a lattice 
interval that is driven by independent birth and death processes at its boundaries 



corresponding to different chemical potentials. Take A„ = {—n, —n + 1, . . . ,n — l,n} 
and 1] G {0, 1}^" a particle configuration evolving with generator 

Gnfiv) = E [/(^^•^+^) - /(^)] 
x=—n 

+ A[e^^''-"(/(r")-/(r^)) + e^^''"(/(r^")-/(r^))] (1.8) 

The first term corresponds to symmetric hopping with exclusion; the two last terms 
are giving birth and death to particles at the ends of the interval with parameters 
hi, h2- One can think here of particle reservoirs, to the left of the system with density 
1/(1+6^^) and to the right with density 1/(1 +6^^). For A = the system is uncoupled 
from the reservoirs and it has all uniform product measures as reversible measures 
with vanishing mean entropy production. For X ^ 0, hi ^ h2 this detailed balance is 
lost and we have positive mean entropy production. Yet, it remains of order unity, 
uniformly in the size n meaning that the mean entropy production density vanishes 
in the thermodynamic limit. This is an instance of a more general fact for interacting 
particle systems that will also be treated in the next section: you cannot by driving 
the system at its boundaries break the time-reversal invariance in the limiting infinite 
volume process, see Proposition 2 below. 

We do not know whether there exists a time-reversal vr for which ( |1.8| ) would give 
rise to generalized detailed balance. 

In this paper we show more generally how breaking of detailed balance is strictly 
equivalent with non-zero mean entropy production. There is no way to get a current 
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and at the same time to have no dissipation (zero mean entropy production). 

In the next section we describe our class of models and we state our main result. 
In section 3, we discuss this result and we give some more background information 
concerning entropy production, reversibility and time-reversal. Section 4 is devoted 
to the proofs. 



2 Models and main result 

2.1 Dynamics 

This subsection describes the assumptions and introduces the necessary notation. 
The configuration space is Q = 5"^ where 5" is a finite set and Z is the regular 
0?— dimensional lattice. Let tt be an involution on O. A special but important case 
is when tt =identity. We assume here that tt commutes with lattice translations 

Let Vo C Z'^ be a finite cube containing the origin and write Vq for any specific non- 
empty set of transformations C/q on ^ satisfying, for every C/q £ 'Po, and for every 

i. (C/oa)(y) = (7(y),foryeyo^ 

ii. C/q-' e 7^0, 

iii. 'kVq'k — Vo, 



iv. If Uo 7^ Uq and Uoa ^ cr, V'^a ^ a then UqG ^ U'^a (for convenience only.) 

We consider the translations = {y -\- x : y & Vq\ and = tJJqT^x to generate 
a dynamics via local translation invariant rates c{Ux-i cr) for the transition a U^cr. 
We assume: 

V. Positivity: c{Uo, a) = when Uoa = a and if not, c{Uo, a) > 0, 

vi. Finite range: there is a finite A G such that for all (7,77 e fl, and Uq G Vq- 

c{Uo,a) = c{Uo,a]^r]xc), 

vii. Translation invariance: for all x e Z'', Ux G T^i, a G fl: c{Ux,<j) = c{Uo,T-x(j) 

The generator L corresponding to the given rates is now defined on local functions / 
as 

That is, a is changed to 77 at rate c(C4, cr) if = UxCr. We will always write p for a 
translation invariant stationary measure for this dynamics. It can be different from 
pir but we assume that also pir is stationary. Finally, p and pir give positive weight 
to all cylinders and writing p^ = Up, we always assume that dp^° / dp{a) > c > 0, 
which, even in the present rather general set-up, can be expected quite generally. 

For Vq — {0} and S — {+1, —1}, the choice UxCr — corresponds to a spinflip 
process. Taking Vq = {0, ei, 62, . . . , e^} with Cq the lattice unit vectors, we can make 
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a spin exchange process or hopping dynamics. We refer to for further details on 
constructing the infinite volume process. 



2.2 Mean Entropy Production 

Put A„ = [—n^nf n 'Sf for large n and define Aj^ as the maximal subset of A„, such 
that for all a; G Aj^ and G P^, ciU^^ cr) depends only on coordinates inside A„, and 
14 C A. Consider now the Markov chain on S*^" with generator 

L„/(a)= E c(t/.,a)[/(f/.a)-/(a)] (2.10) 

and started from a probability measure p_T on S^^ at time -T. The measure at time 
T is denoted by px- Via a Girsanov formula this dynamics gives rise to a Hamiltonian 
(or action functional) on space-time trajectories u (as in |^ 0), with corresponding 
relative energy with respect to time- reversal given by the entropy production (O) 
and here equal to 

i?,(L„,p_T,T)(^) = In p_t(^(-T))- In pT(w(T)) + A5e(^) (2.11) 

with 

+ / [c{U:,,'Ru{s))-c{U^,u{s))]ds (2.12) 

J-T 

where N^'^^co) = J2-T<s<t ^ i^i-'') ~ Uxi^^is")) ^ uj{s~)) is the number of times the 
transformation Ux appeared in the realization u up to time t G [— T, T]. The ex- 
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pression ( |2.12| ) must be interpreted as the variable entropy produced in the reser- 



voirs (environment) when the microscopic system configuration moves from uj{—T) 
to uj{T): To get the total variable entropy production (|2.11|) one should add to ( |2.12D 



the corresponding change in the system's entropy, that are the first two terms in 
( |2.11|) . However, when taking steady state averages, this part vanishes (the entropy 



of the stationary system does not change on average). We can therefore define the 
mean entropy production for the interacting particle system as 

MEP.(L,p) ^ hm hm ^^^e;'^(A5,) (2.13) 



lE^'"^ denotes the expectation with respect to the path space measure, in the station- 
ary distribution p, restricted to trajectories within S^". In other words, the mean 
entropy production is the expectation of the time-reversal breaking part in the space- 
time action functional governing the dynamics. We refer to |^ for a mathematical 
discussion on the existence of the limit (|2.13| ) and for a proof of its non-negativity. 



We refer to 0, ||, Q and Section 3 for further background. 
2.3 Results 

The main question is to see whether for a dynamics where the time-reversal symmetry 
is explictly broken (in the sense that there is no detailed balance), there still can be 
zero mean entropy production (dissipationless steady state). Our main result says 
that this is impossible. 
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Main Theorem: Under the conditions above, MEPt^{L, p) = MEPt^{L, pn) = 
implies that the dynamics satisfies (generalized) detailed balance in the sense that for 
all Uq 

c(nUo\,TTUoa)^^(a) = c(Uo,(t) p - a.s. (2.14) 
dp 



Note that (|2.14| ) is really the analogue of (|1.3|) . Observe also here that ( |2.14| ) implies 



that the densities dp^°/dp are invariant under replacing p by pyr. This follows from 
rewriting (|2.14|) from right to left with a nUocr and Uq ttUq'^it: 

dp^^o^'^ 

c^TTUQ^iijTiUoa) = c^nnUo'K'Kj'K'KUQ^miUoa) — {nUocr) 

ITT -.dpT^i ^ 



and comparing it with the original ( p.l4| ). 



We call Vq complete if every local transformation h : Vt Vt can be written as a 
composition of Ux'- i.e., ii h = U^^ ■ ■ ■ Ux„ for some xi, . . . , x„ G Z'^. 
Corollary 1: If MEP^^lL, p) = MEPt^IL, pir) = and if Vq is complete and vr is 
continuous, then p is a reversible Gibbs measure for the dynamics defined above. 

In the converse to these results was already shown: Suppose that the rates 
satisfy 

c([/,, a) = ci-nU^^-K, TiUxa) exp{-H{Uxa) + H{a)). (2.15) 
This is again the analogue of ( |l.tj| ). The energy difference in ( |2.15| ) should be inter- 
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preted in terms of an absolutely convergent sum of potentials: 

H{aAr]Ac)-H{UVM)= ^ {ViA, ctat^a.) - V{A, Uvm)) , (2.16) 

where {V{A, ■) : 5*^ {—oo, +00), A finite subsets of Z*^), is a translation invariant 
(uniformly) absolutely summable potential: 

^max\V{A,a)\ < +00 (2.17) 

Then, 

MEP^(L,p) = MEP^(L,p7r) = 
When we combine the above we obtain a final 

Corollary 2: Under the conditions of Corollary 1, if there is one translation in- 
variant stationary measure p for which p = pir and MEP^ (L,p) = 0, then also 
MEPt^IL,!/) = for all translation invariant stationary measures v and they are all 
Gibbsian for the same potential. 

A caveat in the above main result is to understand better the relation between 
MEP^(L,p) and MEP^(L,p7r). To this we can only add that MEP^(7rL7r, p) = 
MEP7r(L, p7r), as can be verified from a direct computation starting with ( [4.24| ). 



The simplest illustration of all this was already obtained in ^ for a spinfiip process. 
Example B, (|1.7|) , deals with a spinfiip process but there the time-reversal vr does 
not commute with translations. As will be seen from the proof, that is indeed not 
essential as long as the dynamics and the stationary measure are translation invari- 
ant. Of course, one should then be extra careful with condition iii. but also this can 
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be modified accordingly. It will also be clear that more general lattice structures and 
configuration spaces can be employed (e.g. already in Example A). 
Finally, for completeness we come back to the situation of Example C in Section 
1.1. For this we must leave the translation invariant infinite volume context and ask 
whether boundary driven interacting particle systems can give rise to non-vanishing 
mean entropy production density in the thermodynamic limit. The question can be 
formalized as follows. We consider a process on S^" with generator Gn generalizing 



ryes-* 

^ C A„, \ Al 
diamA < r 

where o"^''' = (Ta'^t]a equals a outside the set A which has a diameter (maximal lattice 
distance within) less than a given constant r. 

Here the generator L„ is given by ( p.lO| ) but with rates verifying condition ( p.l5| ) for a 



finite range potential, and rates A;^^(cr, 77) as in ( |1.8D inducing configurational changes 
at the boundary of A„. We further assume that the /c^"* (a, 77) are uniformly bounded 
from below and from above. In other words, we have a bulk dynamics generated 
by Ln with rates satisfying (generalized) detailed balance, and at the boundary the 
configuration can change quite arbitrarily (but in a local and bounded way). We 
suppose that p„ is the unique stationary measure of this dynamics and for simplicity 
we only treat the case tt = id. We are interested in the mean entropy production 

15 



MEP{Gn,pn) defined in (^) (with n = id). 

Proposition 2: There is a constant K so that MEP^(G'„,p„) < Kn''-~^ 
The proofs of the above results are postponed to Section 4. 

3 Discussion 

We briefly discuss some concepts that are important for our result. 
3.1 Time- reversal 

By this we usually mean a transformation on phase space VL which, for a many-particle 
system, is defined particle-wise or, for spatially extended systems, is sufficiently lo- 
cal. Physically speaking, its precise nature follows from kinematical considerations 
on the dynamical variables. In classical mechanics, it reverses the momenta of all 
the particles but in the presence of say an electromagnetic potential, considered part 
of the system, one can add an extra transformation reversing also the magnetic field 
and thus making the Lorentz force time-reversal invariant. In our case, we have a 
configuration space Vt = S with Z the (i-dimensional lattice and S a finite set. 
Time-reversal is an involution vr on fi, tt^ = id. Time-reversal extends to a transfor- 
mation on path-space, as introduced for (|1.4|) , by reversing the trajectories. That 
is, if we have a trajectory {uit^t G [— T, T]) then the time-reversed trajectory 6'^(a;) is 
given by (61^ (a;))* = nuj-f 
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3.2 Reversibility 

Dynamic reversibility is a property of the dynamics itself under time-reversal. It says 
that if one trajectory u) of the system is possible, so is its time-reversed Qt,{uS). For a 
deterministic system where cut = 0(it)ci;o with an invertible flow on phase space, 
it says that 4>{t)~^ = 7r0(t)7r, that is a symmetry that anticommutes with the time 
evolution. For a stochastic dynamics this is implied by assuming that if a transition 
a ^ Ua is possible (positive transition rate), then also the same is true for its time- 
reversal vrf/cr — >• TTO". 

Microscopic reversibility is a consequence of dynamic reversibility in case of an equi- 
librium dynamics. For our purposes here we do not make a distinction with the 
condition of detailed balance. When the dynamics is driven away from equilibrium, 
the resulting stochastic model will not satisfy detailed balance. Usually this produces 
a current in the system (but that need not be true in general, see an example in 0). 
On the other hand, a net current signifies the breaking of the detailed balance condi- 
tion. In general we like to distinguish between two classes of finite volume dynamics 
where microscopic reversibility is explicitly broken. These are boundary driven ver- 
sus bulk driven dynamics depending on the extensivity of the perturbation from an 
equilibrium dynamics. In the bulk driven case, one usually verifies so called local 
detailed balance, i.e., (|2.15|) is changed into 



c(f/,, a) = c{nU-\, nU^a) exp{-H{U^a) + if((T))e^*(^^'^''^) 
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where E is some amplitude of an external field and $ is antisymmetric, $(7r?7, ira) = 
— $(cr, 1]), see e.g. [0. Note also that then, necessarily, the relative energies HiUxc) — 
H[a) are invariant under exchanging H with Ht^. 

In boundary driven systems, the process becomes non-translation invariant and the 
rates remain of the form (|2.15|) in the bulk (that is for x well inside the considered 



finite volume) while more or less arbitrary on the boundary. This was the case for 
Example C in Section 1.1 and was formalized for Proposition 2. Note that there 
is in fact an example of a boundary driven system where uniformly in the size of 
the system a bulk current can be maintained. This is the nonequilibrium harmonic 



crystal treated in [^, Q where the heat flux is proportional to the boundary temper- 
ature difference rather than to the temperature gradient (infinite heat conductivity 
in the thermodynamic limit). Such 'superconductors' do not exist in the context of 
interacting particle systems as discussed in the present paper. 



3.3 Entropy production 

In phenomenological thermodynamics, entropy production appears in open driven 
systems as the product of thermodynamic fluxes and forces. The forces are gradients 
of intensive quantities (like chemical potential) generating the currents. The entropy 
production is identified from a balance equation for the time-derivative of an entropy 
density which is defined in systems close to equilibrium. The definition of entropy 
production as we use it here in statistical mechanics comes from p, |, H, |^, [11], |12|, [14 
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and we refer to the review The mean entropy production appears there and in 
(|1.4|)-( pT5D as a relative entropy (density) for the process with respect to its time- 
reversaL That immediately invites the following thought (we are grateful to Senya 
Shlosman for pointing to this): In equilibrium statistical mechanics, if two translation 
invariant Gibbs measures have zero relative entropy density, then they must both be 
Gibbsian for the same interaction potential (but not necessarily equal e.g. because 
of spontaneous symmetry breaking). Apply this to the space-time measures obtained 
for the process Pp and the time-reversed process Pp0 as introduced for (p^). Here 
we take vr = identity to avoid extra complications. In some sense, both processes are 
Gibbs measures. Thus, if the mean entropy production is zero, then the process itself 
and its time-reversal have the same (space-time) action functional. Because they 
also have the same marginals p, they must in fact coincide (hence no spontaneous 
time-reversal breaking). Hence, zero mean entropy production implies microscopic 
reversibility. While convincing on a superficial level, unfortunately the details of this 
argument are technically cumbersome and a direct sufficiently general proof along 
this line has not been found. 

The only more recent paper that we know of concerning time-reversal symmetry 
and the relation with entropy production is [|l|. The set-up there is however quite 
different from ours. Time-reversal symmetry is there associated with the anticom- 
mutation of an involution with the time evolution, what we have called dynamic 
reversibility in the above. In our discussions here, we deal with spatially extended 
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stochastic dynamics and the breaking of microscopic reversibihty. 



4 Proofs 



Lemma 1: Under the conditions of Section 2.1, for a translation invariant stationary 



measure u. 



J du{a)c{Uo,a)log-^{a) = Q 



(4.18) 



Proof: Let J^a be the a- field generated by (T^,x G A. Denote by u\, respectively 
the restrictions of v and . Then we have 



dUA 



diy^° 
dv 



Ta 



Since dv^° /dv e L^{du) for all C/q, we find using the martingale convergence theorem 
that 



lim ^ - ^ 
j^^jf dvx dv 



(4.19) 



in L^{dv). Let z/ be the product measure on Q having as marginals the uniform 
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measure on S. Prom stationarity applied to the local function f\ — du^/ di>\ we find 



di/A 

du^'' 
dv 



dUA 



du{a)c{U,, a)[\og^ -log 
= ^ ^ / diy{a)c{U^, a) log 
= ^ / diy{a)c{U^, a) log 

+ J2 I dv{a)c{U:,,a) 



= |A'| ^ / dK<7)c(C/o,<7)log— ^(<7) 

f/oe-Po ^ 



log log — — 



The last equality uses translation invariance. We have used the notation A' = {x e 
n A 7^ 0} and the expression 



dv 



We thus have 



5^ fdi.{a)c{Uo,a)log^{a)\ < E E I f di^{a)c{U^,a)Fj^^{a)\ 

^ ^mE E / ^Hi^r°J, (4.20) 



a;GA' [/oSPo ' 

by the translation invariance of and M bounds the rates. Now we use the general 
fact that if converges to / in L}{dv) and both /„, / are bounded from below by 
some constant c > 0, then log/„ converges to log / in L^{dv). This fact implies that 
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for any given e > 0, we can choose A C such that for all A' D A: 



max/dHfi;?|<^, with m^N 



Choose now A C Z'^ so large that 

|{x G A + xn ^ 0}| 

W\ " 2MNsv.Vw,uA\Fw'\Widu) 

We then conclude that 



< 



fduia) c{Uo,a) log^(a)|<^ 



,,,, |{x G A' : A + xnA^ ^ 0}| 
MN ^ sup I I 

|A I w,Uo 



< e (4.21) 



Proof of Main Theorem: Define 



c{Uo, a) = c{tt, p; Uq, a) = c{nUo\, ^^oo-)-^(o-) 



and substitute it in 



We get four terms, ( [4.22|) 



J2 [ dpia) [c(t/o, a) - c(t/o, a)] log (4.22) 



J2lfdpi-)^og-^l^l^\— + [ dp{a)c{U„ a) log 



+ / dp^°(a)c(7r[/o-V, vr^oa) log ''.^"''^ + / cip^°(a)c(7rt/o-V, vrt/oa) log ^ 
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The second term is zero by Lemma 1. The fourth term is also zero because, using 
condition iii, we can change txUq^h Uq in the sum over Vq getting it equal to 



^ J dpTT{a)c{Uo, a) log 



dipir) 



Uo 



which is zero, again by Lemma 1 applied to the stationary measure pn. Again usin^ 
iii, we can also rewrite the third term as 



E fdpA^) log , ^(^"'^^ , 



UoeVo ' 

Therefore, what remains of ( [1.23| ) is the sum of the first and the third term so that 
equals 



We now recall that the mean entropy production ( p.l3| ) equals 



c{Uo,cr) 



c{ttUq TT, TT f/oO") 



+ / rfp(a)[c([/o,7ra)-c([/o,a)]) (4.24) 



This was derived from ( ^.12| ) in We conclude therefore that ( |4.22| ) equals MFjPt^(L, p) + 



MEP^(L, pre) which is zero by hypothesis. This implies the statement of the Theorem.i 



Proof of Corollary 1 Since the Radon-Nikodym derivative of p ° with respect to 
p is a local function for all Uq and since by assumption, we can generate with the Uq 
all local excitations a' from a, it means that p has a continuous version for its local 
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conditional distributions. 



Proof of Corollary 2 From the main result and Corollary 1 it follows that p is 
a translation invariant stationary Gibbs measure and (|2.15|) must be satisfied. All 



other translation invariant stationary measures must be Gibbsian and for the same 
potential, see e.g. 0. From the results in as cited above the statement of Corol- 
lary 2, it follows that every other stationary translation invariant measure must have 
zero mean entropy production. ■ 



Proof of Proposition 2 From the definition ( |1.5| ) we must first compute the relative 



action under time reversal, that is 



This can be done via a Girsanov formula and we obtain the analogue of (|2.12|) : 



A,ij,n I 



diam(A) < r 

The first integral is really a sum over all the times when the trajectory makes a jump 
from the action of one of the Ux'-, the second integral is a sum over all times when a 
configuration a is replacing 1^(5^) in a set A on the boundary. In order to further 
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clarify this formula, let us first look at trajectories where no boundary transitions 
take place (or, what amounts to the same, take /c = for the moment). Then, we 
only keep the first term, that is just (|2.12|) , in case tt = id: 

y y flog ^(^-^(^")) dN^. 



-T 



But if we insert the detailed balance condition (|2.15|) , the above expression telescopes 



to 

H{uj{-T))-H{uj{T)) 
and the mean entropy production is zero by stationarity. 

Turning to the general case we let {s,}^^^ be the set of times at which boundary 
transitions occur in the sets Aj, i = 1, . . . , g, for the trajectory uj. These are random 
but we fix them as — T < si < S2 < .. < Sg < T. The important thing to realize now 
is that while the perfect telescoping of above is broken at each of these times, it can 
be restored by adding and subtracting. More precisely, we have 

Rn{uj) = H{uj{~T))~H{uj{si)) + H{uj{s,))-H{uj{s^)) + ... 
+ HiLuis,))-HiiuiT)) + \og f'^ + 

+ log -r4 ^ h . . . + log — 

But by the absolute convergence of the interaction potential we have 

\H{uj{si)-H{uj{s,))\<rC 
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for some constant C, since uj{s~) and uj{si) only differ in the set Aj. Therefore the 
telescoping of the terms involving energy differences can be restored upon inserting 
q terms of order unity. 

As for the other terms, we have assumed uniform boundedness so that we get 

M 

|i4HI <g(rC + log— ) 

where M and e are constant upper and lower bounds for the transition rates /c*^"-* . As 
the expectation oi q — q{u!) under lE^^ is proportional to Tl^A^I, the proposition is 
proved. ■ 

Acknowledgment: CM. thanks Senya Shlosman for some very useful discussions. 
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